MULTIVARIATE FUSS-CATALAN NUMBERS 

J.-C. AVAL 

Ji—^ I Abstract. Catalan numbers C{n) = ;nr[(^) enumerate binary trees and Dyck 

f~^ ' paths. The distribution of paths with respect to their number k of factors is given 

CSI . by baUot numbers B{n,k) = ^ij^("^ )■ These integers are known to satisfy sim- 

^ ' pie recurrence, which may be visualised in a "Catalan triangle" , a lower-triangular 

Q , two-dimensional array. It is surprising that the extension of this construction to 

1/'^ ' 3 dimensions generates integers B^{n,k,l) that give a 2-parameter distribution of 

Csin) = 2n+i ('^n ) ' 'which may be called order-3 Fuss-Catalan numbers, and enu- 
merate ternary trees. The aim of this paper is a study of these integers B^{n, k, I). 
We obtain an explicit formula and a description in terms of trees and paths. Finally, 
^^ ' we extend our construction to p-dimensional arrays, and in this case we obtain a 

r) . (p — l)-parameter distribution of Cp{n) ~ 7— zjiti+t ('r") ' ^^^'^ number of p-ary trees. 

^" 



\o 



1. Catalan triangle, binary trees, and Dyck paths 

We recall in this section well-known results about Catalan numbers and ballot 
\0 [ numbers. 

5 ' The Catalan numbers 

C[n) 



> 



0^ 
O 



r- 



X 



c^ 



n + I \ n 



are integers that appear in many combinatorial problems. These numbers first ap- 
O ■ peared in Euler's work as the number of triangulations of a polygon by mean of 

non-intersecting diagonals. Stanley [T6l [T7j maintains a dynamic list of exercises re- 
lated to Catalan numbers, including (at this date) 127 combinatorial interpretations. 



^ ■ Closely related to Catalan numbers are ballot numbers. Their name is due to 



the fact that they are the solution of the so-called ballot problem: we consider an 
election between two candidates A and B, which respectively receive a > b votes. 
The question is: what is the probability that during the counting of votes, A stays 
ahead of B? The answer will be given below, and we refer to [5] for Bertrand's first 
solution, and to [1] for Andre's beautiful solution using the "reflection principle" . 

Since our goal here is different, we shall neither define ballot numbers by the previ- 
ous statement, nor by their explicit formula, but we introduce integers B{n, k) defined 
for a positive integer n and a nonnegative integer k by the following conditions: 

• 5(1,0) = 1; 

• Vn > 1 and < A; < n, B{n, k) = X;J=o ^(^ -1^0; 

• Wk>n, B{n, k) = 0. 
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Observe that the recursive formula in the second condition is equivalent to: 



:i-i) 



B{n, k) = B{n - 1, A;) + B{n, k-l] 



We shall present the B{n,k)^s by the following triangular representation (zero 
entries are omitted) where moving down increases n and moving right increases k. 
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1 
1 
1 
1 
1 



1 
2 
3 

4 
5 



2 

5 5 

9 14 14 

14 28 42 42 



The crucial observation is that computing the horizontal sums of these integers give 
: 1, 2, 5, 14, 42, 132. We recognize the first terms of the Catalan series, and this 
intuition will be settled in Proposition 11.11 



after introducing combinatorial objects. 

A binary tree is a tree in which every internal node has exactly 2 sons. The 
number of binary trees with n internal nodes is given by the n-th Catalan number. 
The nodes of the following tree are labelled to explain a bijection described in the 
next paragraph: internal nodes are labelled by letters and external nodes by numbers. 




A Dyck path is a path consisting of steps (1, 1) and (1, —1), starting from (0,0), 
ending at (2n, 0), and remaining in the half-plane y > (we shall sometimes say 
"remaining above the horizontal axis" in the same sense). The number of Dyck 
paths of length 2n is also given by the n-th Catalan number. More precisely, the 
depth-first search of the tree gives a bijection between binary trees and Dyck paths: 
we associate to each external node (except the left-most one) a (1, 1) step and to each 
internal node a (1,-1) step by searching recursively the left son, then the right son, 
then the root. As an example, we show below the Dyck path corresponding to the 
binary tree given above. The labels on steps correspond to those on the nodes of the 
tree. 

A A ^ 




An important parameter in our study will be the length of the right-most sequence 
of (1, —1) of the path. This parameter equals 2 in our example. Observe that under 
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the correspondence between paths and trees, this parameter corresponds to the length 
of the right-most string of right sons in the tree. We shall use the expressions last 
down sequence and last right string, for these parts of the path and of the tree. 

Now we come to the announced result. It is well-known and simple, but is the 
starting point of our work. 

Proposition 1.1. For n a positive integer, we have the following equality: 

n-l 






fc=0 

Proof. Let us denote by Cn,k the set of Dyck paths of length 2n with a last down 
sequence of length equal to n — k. 

We shall prove that -B(n, k) is the cardinality of Cn,k- 

The proof is done recursively on n. If n = 1, this is trivial. If n > 1, let us suppose 
that B{n — 1, k) is the cardinality of Cn-i,k for < /c < n — 1. Let us consider an 
element of Cn,k- If we erase the last step (1,1) and the following step (1,-1), we 
obtain a Dyck path of length 2(n — 1), with a last decreasing sequence of length 
n — 1 — I with I <k. If we keep track of the integer fc, we obtain a bijection between 
Cn k and UKfcC„,-_i I. We mention that this process is very similar to the ECO method 

This is a combinatorial proof of Proposition 11.11 D 

Remark 1.2. The integers B{n,k) are known as ballot numbers and are given by 
the explicit formula: 

, , , a — h fa + h 
(1.2) B{a,h) 



a + b 

This expression can be obtained shortly by checking the recurrence (ILip . We can 
alternatively use the reflection principle (see [11] for a clear presentation), or the 
cycle lemma (c/. [6]), which will be used in the next sections to obtain a formula in 
the general case. 

The expression (11.21) constitutes a solution to the ballot problem. For this, we 
use a classical interpretation in terms of paths: we represent a vote for A by an up 
step, and a vote for B by a down step. The total number of countings of votes, or of 
paths from (0, 0) to {a + b, a — b) , is given by the binomial ('^^ ) . The countings such 
that A stays ahead of B correspond to paths remaining above the horizontal axis. 
Their number is given by B{a,b). This implies that the probability asked for at the 
beginning of this section is ^^. 

Of course , we could have used expression (11.21) to prove Proposition 1 1 . 1 1 by a simple 
computation, but our proof explains more about the combinatorial objects and can 
be adapted to ternary trees in the next section. 

It should be mentionned here that our study of multivariate Fuss-Catalan has 
nothing to do with the many-candidate ballot problem, as considered for example 
in [18] or [13]. In particular, the multivariate ballot numbers considered in these 
papers do not sum to the Fuss-Catalan numbers Cp{n) = , _iL_|_i (^^)- Conversely, 
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the numbers studied in the present article do not give any answer to the generahzed 
ballot problem. 

Remark 1.3. Our Catalan array presents similarities with Riordan arrays, but it is 
not a Riordan array. It may be useful to explicate this point. We recall (c/. [H]) that 
a Riordan array M = {rriij) G C^^^ is defined with respect to 2 generating functions 

fi'(^) = X^ dnx"" and /(x) = ^ /„a;" 

and is such that 

Mj{x) = ^mn,jx'^ = g{x)f{xy. 

n>0 

We easily observe that a Riordan array with the first two columns Mq{x) and Mi{x) 
of our Catalan array is relative to g{x) = j^ and f{x) = j^, which gives the Pascal 
triangle. 

In fact, the Riordan array relative to g{x) = C(x) = ^ C{n)x"' and f{x) = x C{x) 
gives the ballot numbers, but requires knowledge of the Catalan numbers. 

2. Fuss- Catalan tetrahedron and ternary trees 

2.1. Definitions. This section, which is the heart of this work, is the study of a 
3-dimensional analogue of the Catalan triangle of the previous section. That is we 
consider exactly the same recurrence, and let the array grow, not in 2, but in 3 di- 
mensions. More precisely, we introduce the sequence B^i^n, k, I) indexed by a positive 
integer n and nonnegative integers k and I, and defined recursively by: 

. 53(1,0,0) = 1; 

• yn> 1, k + l <n, B3{n,k,l) = Y.o<i<k,o<j<i ^sin - l,i, j); 

• \/k + l>n,B3{n,k,l)=0. 

Observe that the recursive formula in the second condition is equivalent to: 

(2.1) 53(n, k, I) = B^in - 1, k, l) + Bs{n, k-1, l) + B3{n, k, I -I) - B^in^k-lJ-l) 

and this expression can be used to make some computations lighter, but the presen- 
tation above explains more about the generalization of the definition of the ballot 
numbers B{n,k). 

Because of the planar structure of the sheet of paper, we are forced to present the 
tetrahedron of B^^n, k, /)'s by its sections with a given n. 



n 



n 



n 



n 
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n 



1 


4 


9 


14 


4 


15 


30 


35 


9 


30 


45 




14 


35 






14 









14 



It is clear that Bs{n,k,0) = B3{n,0,k) = B{n,k). The reader may easily check 
that when we compute ^^; _B3(n, /c, /), we obtain: 1, 3, 12, 55, 273. These integers 
are the first terms of the following sequence (c/. [T5]): 

1 /3n' 
Cain) 



2n + 1 \n 
This fact will be proven in Proposition 12. 1[ 

2.2. Combinatorialinterpretation. Fusqj-Catalan numbers (c/. [9l[TT]) are given 
by the formula 

1 / pn' 



(2.2) 



Cpin) 



{p — l)n + 1 



n 



and C^{n) appear as order-3 Fuss-Catalan numbers. The integers 6*3(71) are known 
to count ternary trees, ie. trees in which every internal node has exactly 3 sons. 




Ternary trees are in bijection with 2-Dyck paths, which are defined as paths from 
(0,0) to (3n, 0) with steps (1,1) and (1,-2), and remaining above the line y = 0. 
The bijection between these objects is the same as in the case of binary trees, ie. a 
depth-first search, with the difference that here an internal node is translated into 
a (1, —2) step. To illustrate this bijection, we give the path corresponding to the 
previous example of ternary tree: 




We shall consider these paths with respect to the position of their down steps. 
The height of a down step is defined as the height of its end-point. Let T>n,k,i denote 
the set of 2-Dyck paths of length 3n, with k down steps at even height and / down 



Nikolai Fuss (Basel, 1755 - St Petersburg. 1826) helped Eulcr prepare over 250 articles for 
publication over a period on about seven years in which he acted as Euler's assistant, and was from 
1800 to 1826 permanent secretary to the St Petersburg Academy. 
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steps at odd height, excluding the last seqence of down steps. By definition, the last 
sequence of down steps is of length n — k — I. 



Proposition 2.1. We have 



Vn>0, ^B3{n,k,l) = C3{n) 

k,l 



2n + l\ n 



3n 



Moreover, B^in^k^l) is the cardinality ofV^. 



k,l- 



Proof. Let k and / be fixed. Let us consider an element of Vn^k,!- If we cut this path 
after its {2n — 2)-th up step, and complete with down steps, we obtain a 2-Dyck path 
of length 3(n — 1) (see figure below). It is clear that this path is an element of ^n,ij 
for some i < k and j < I. We can furthermore reconstruct the original path from the 
truncated one, if we know k and /. We only have to delete the last sequence of down 
steps (here the dashed line), to draw k — i down steps, one up step, / — j down steps, 
one up step, and to complete with down steps. This gives a bijection from Vn^k,! to 
^o<i<k,o<j<i1^n-i,i,j, which implies Proposition [2111 




D 



Remark 2.2. It is interesting to translate the bi-statistic introduced on 2-Dyck 
paths to the case of ternary trees. As previously, we consider the depth-first search 
of the tree, and shall not consider the last right string. We define 7^,fc,z as the set of 
ternary trees with n internal nodes, k of them being encountered in the search after 
an even number of leaves and / after and odd number of leaves. By the bijection 
between trees and paths, and Proposition 12.11 we have that the cardinality of Tn^k,i 
is -83(72, k, I). 
Remark 2.3. It is clear from the definition that: 

Bs{n,k,l) = B3{n,l,k). 

But this fact is not obvious when considering trees or paths, since the statistics 
defined are not clearly symmetric. To explain this, we can introduce an involution on 
the set of ternary trees which sends an element of Tn^k,i to 7^,z,fc- To do this, we can 
exchange for each node of the last right string its left and its middle son, as in the 
following picture. Since the number of leaves of a ternary tree is odd, every "even" 
node becomes an odd one, and conversely. 
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2.3. Explicit formula. Now a natural question is to obtain explicit formulas for the 
i?3(n, k, I). The answer is given by the following proposition. 

Proposition 2.2. The integers B^{n,k^l) are given by 



(2.3) 



B3{n,k,l) 



n + k 
k 



n + l 
I 



n 



k-l 



n 



k 



Proof. We use a combinatorial method to enumerate Vn^kj- The method is a variation 
of the cycle lemma [6] (called "penetrating analysis" in |12j). 

If we forget the condition of "positivity" {ie. the path remains above the line 
y = 0), and cut the last down sequence, a path consists in: 

• at even height: n up steps, and k down steps; 

• at odd height: n up steps and / down steps. 

An important remark is to remember that an element of 'Dn,k,i has an up step just 
before the last sequence of down steps! If we suppose that all these steps are distin- 
guished, we obtain: 



(n+fc)! 

nlkl 
(n-l+Q! 



choices for even places; 



- (j^_iYi, choices for odd places (we cannot put any odd down step after the 
last odd step). 

Now we group the paths which are "even permutations" of a given path P. By even 
permutations, we mean cycle permutations which preserve the parity of the height 
of the steps. 

We want to prove that the proportion of elements of T>n^k,i in any even orbit {i.e. 
in any orbit under even permutations) is given by ""^f ~' . 

We suppose first that the path P is acyclic (as a word): P cannot be written as 
P = U^ with U a word in the two letters (1, 1) and (1, —2). It is clear that such a 
path gives n + k different even permutations. Now we have to keep only those which 
give elements of T>n^k,i- To do this we consider the concatenation of P and P', which 
is a duplicate of P. 
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The cyclic permutations of P (not necessarily even) are the subpaths oi P + P' 
of horizontal length 2n + k + I. The number of such paths that remain above the 
horizontal axis, and end with an up step, is the number of (up) steps of P in the 
light of an horizontal light source coming from the right. The only transformation is 
to put the illuminated up step at the end of the path. The number of illuminated up 
steps is 2n — 2k — 21, since every down step puts two up steps in the shadow. Among 
these 2n — 2k — 21 permutations, only half are even (observe that the set of heights of 
the illuminated steps is an interval). Thus n — k — I paths among the n + k elements 
of the orbit are in 'Dn,k,i- 

Now we observe that if P is p-cyclic, then its orbit has p times less elements, and we 
obtain p times fewer different paths, whence the proportion of elements of elements 
of Pn.fc,/ in this orbit is: 

{n — k — l)/p n — k — I 
(n + k)/p n + k 

Finally, we obtain that the cardinality of Pn,fc,/ is 

{n + k)\{n-l + iy.n-k-l 
n\k\ (n-1)!/! n + k 

which was to be proved. 

D 

Remark 2.5. The equation (12. 3p is of course symmetric in k and /: 

^n + k\ fn + I — l\n — k — I fn + k — l\ fn + 1 — \\n — k — I 



k J \ I J n + k \ k J \ I J n 

Remark 2.6. I made the choice to present a combinatorial proof of Proposition 12.21 
The interest is to show how these formulas are obtained, and to allow easy gener- 
alizations (c/. the next section). It is also possible to check directly the recurrence 

(EU). 

2.4. Generating function. Let F denote the following generating series: 
(2.4) F{t,x,y) = l+ Yl B,{n,k,l)ex''y'. 

0<k+l<n 

In this expression, the term "1" corresponds to the empty 2-Dyck path. Thus F is 
the generating series of 2-Dyck paths with respect to their length (variable t), their 
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number of down steps at even height -excluding the last down sequence- (variable 
x), and their number of down steps at odd height (variable y). 
We also introduce 

(2.5) G{t,x,y) = l+ Yl 53(n,fc,/)rx"-V, 

0<k+l<n 

ie. G is the generating function of 2-Dyck paths with respect to their length, to the 
number of down steps at even height -including the last down sequence- and to the 
number of down steps at odd height. 

To obtain equations for F and G, we decompose a non-empty 2-Dyck path by 
looking at two points: a, defined as the last return to the axis (except the final point 
(3n, 0)), and (3 defined as the last point at height 1 after a. This gives the following 
decomposition of a 2-Dyck path 

P = Pi (up) P2 (up) P3 (down), 

with Pi, P2, P3 any 2-Dyck path (maybe empty). 




By observing that the up steps after a and P change the parity of the height, this 
gives the two following equations : 

(2.6) F{t, x,y) = l + G{t, x, y) x G(t, y, x) x t.F(t, x, y) 
and in the same way 

(2.7) G{t, X, y) = 1 + G{t, x, y) x G{t, y, x) x t.G{t, x, y).x. 
By permuting the variables x and y in (12.71) . we obtain 

(2.8) Git, y, x) = 1 + Git, y, x) x G(t, x, y) x t.G(t, y, x).y 

and we can eliminate G{t, y, x) from (12. 7p and (12. Sp to obtain the following result. 
Proposition 2.3. The generating function F of the B^ 's is given by: 

^^^'""'^^ = l-tG{t,x,y)Git,y,x) 
where Git, x, y) is a solution of the algebraic equation 

tx^G^ + iy- x)G'^ + ix- 2y)G + y = 0. 

An alternate approach to the generating function is to use formula 12.31 to obtain 
what MacMahon called a "redundant generating function " (c/. [10]), since it contains 
terms other than those which are combinatorially significant. 
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To do this we extend the recursive definition of B^[n,k,l) as follows: we define 
B'^{n, k, I) for integers n > 0, and k,l > by: 

^n + k — l\/'n + l — l\n — k — I 
k 



(2.9) 



B',in,k,l) 



I I n 

Of course when k + I > n, the integer B'^{n, k, I) is negative. 

As an example, here is the "section" of the array of B'^{n, k, I) with n = 3. 

12 2 0-5 

2 3 0-10 -30 



2 





-12 


-40 


-90 





-10 


-40 


-100 


-200 


-5 


-30 


-90 


-200 


-375 



where c{n, k, I) 



The equation (12.91) is equivalent to the following recursive definition: Vn < or 
kj <0, B'^iji, kj) = and 

B'^{n, k, I) = B'^in-l, k, l)+B'^{n, k-l, l)+B'^{n, k, l-l)-B'^{n, k-l, /-l)+c(n, k, I) 

' +1 if(n,A;,/) = (1,0,0) 

-1 if(n,A;,0 = (0,1,0) 

-1 if(n,A;,/) = (0,0,1) 

+2 if(n,A;,/) = (0,1,1) 

otherwise. 

The interest of this presentation is to give a simple (rational!) generating series. 
Indeed, it is quite simple to deduce from the recursive definition of -83(77., /c, /) that: 

(2.10) ^-- — -^ ' t-x-y^2xy 



J2B3in,kJ)ex'y' 

n,k,l 



1 — t — X — y + xy 
This formula then encodes the B^{n^ k, I) since yk + l < n, B^ln, k, I) = B'^{n, k, I). 

3. Fuss- Catalan ^-simplex and ^-ary trees 

The aim of this final section is to present an extension of the results of Section 
2 to p-dimensional recursive sequences. In the same spirit, we define the sequence 
Bp{n, ki, k2, . . . , kp_i) by the recurrence: 

. 5j,(l,0,0,...,0) = l; 

• Vn > 1 and < ki + k2 + ■ ■ ■ + kp^i < n, 

.,kp_i)= ^ B{n-l,ii,i2,---,ip-i)] 

0<ii<fei,...,0<«p— i<fcp-i 



Bp{n,ki,k2,. ..,k 



• \/ki + k2-\ h kp^i > n, BJn, ki, k2,...,k 



v-i) 



0. 



Every result of Section 2 extends to general p. We shall only give the main results, 
since the proofs are straightforward generalizations of the proofs in the previous 
section. 
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Proposition 3.1. 

The integers Cp{n) are order-p Fuss-Catalan numbers and enumerate p-ary trees, 
or alternatively p-Dyck paths (the down steps are (1, —p))- In this general case, the 
recursive definition of Bp{n, ki,k2, . . . , kp_i) gives rise to p — 1 statistics on trees and 
paths analogous to those defined in section 3. 

Remark 3.2. By the same method as in the previous section, it is possible to obtain 
an explicit formula for these multivariate Fuss-Catalan numbers: 

Bp{n, ki, k2, . . . , kp_i) - ' '^ ' ' '' *=^ * 




v«=l 



Comment. The numbers B^{n,k,l) first arose in a question of algebraic combi- 
natorics (c/. [2]). Let X„ be the ideal generated by 5-quasisymmetric polynomials in 
the 2n variables Xi, . . . , x„ and t/i, ...,?/„ (c/. [1]) without constant term. We denote 
by Q„, the quotient Q[xi, . . . , x„, yi, . . . , yn]/In and by Q^'' the bihomogeneous com- 
ponent of Q„ of degree k in xi, . . . ,Xn and degree / in ?/i, . . . , ?/„. It is proven in [2] 
that: 

1- r>.ki n/ 1 j\ fn + k-l\fn + l-l\n-k-l 
dimQ^'' = 53(n,A;,/) = 



/ / n 

Acknowledgement. The author is very grateful to referees (and in particular 
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